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ABSTRACT 

We test the asymmetry of the Cosmic Microwave Background anisotropy jointly 
in temperature and polarization. We study the hemispherical asymmetry, previously 
found only in the temperature held, with respect to the axis identified by Hansen 
et al. (2009). To this extent, we make use of the low resolution WMAP 5 year tem- 
perature and pola rization iV side = 16 ma ps and the optimal angular power spectrum 
estimator BolPol (|Gruppuso et al.ll2009D . We consider two simple estimators for the 
power asymmetry and we compare our findings with Monte Carlo simulations which 
take into account the full noise covariance matrix. We confirm an excess of power 
in temperature angular power spectrum in the Southern hemisphere at a significant 
level, between 3u and 4<r depending on the exact range of multipoles considered. We 
do not find significant power asymmetry in the gradient (curl) component EE (BB) 
of polarized angular spectra. Also cross-correlation power spectra, i.e. TE, TB, EB, 
sh ow no signific a nt hem ispherical asymmetry. We also show that the Cold Spot found 
by IVielva et al.l ((2004) in the Southern Galactic hemisphere does not alter the sig- 
nificance of the hemispherical asymmetries on multipoles which can be probed by 
maps at resolution N s id. e = 16. Although the origin of the hemispherical asymmetry 
in temperature remains unclear, the study of the polarization patter could add useful 
information on its explanation. We therefore forecast by Monte Carlo the Planck 
capabilities in probing polarization asymmetries. 
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- methods: statistical - cosmology: observations 



1 INTRODUCTION 

Great attention has been devoted to a hemispherical power 
asymmetry in the intensity pattern of the Cos mic Microwave 
Background (CMB) as seen by WMAP l|Hinshaw et all 
l2009l : lDunklev et al.ll2009| ) . Such asymmetry has been origi- 
nally found in WMAP 1st year release and ap pears to lay on 
an ax i s nearly orthogonal t o the ecliptic plane (jEriksen et al.l 
120041 : lHansen et aiTl2004l) . It has been confirmed in the 
WMAP three year and five year release (lEriksen et al 1 l2007l : 
lHansen et all 120091 ; iHoftuft et alj 120091 ) and it is present 
in the COBE data as well, although with lower signifi- 
cance. The temperature power spectra of the opposing hemi- 
spheres are inconsistent at 3a to 4cr depending on the 
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range of multipoles considered. The asymmetry ha s been 
detected in low resolution maps l|Eriksen et al.l l2004). both 
in angular and multipoles space, but it extends to much 
smaller angular sca les in the multipole range 81 — [2, 600] 
l|Hansen et al.l 12009). It is unclear whether this hemispher- 
ical asymmetry is primordial or due to unknown residual 
foreground /systematics . 

Whereas several groups have performed different and in- 
dependent investigations on the CMB temperature pattern, 
the joint analysis of the CMB temperature and polarization 
pattern has not been performed yet. The information en- 
coded in the polarization pattern may turn out extremely 
useful to clarify the presence of the hemispherical asymme- 
try shedding light on its origin. Low resolution WMAP 5 
year maps in (T, Q, U) with relative noise covariance matri- 
ces are publicly available: these public maps have allowed 
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a re-analysis by a quadratic maximum likelihood (hence- 
forth QML) estimator of the low multipol e angular power 
spect rum in temperature and polarization (|Gruppuso et al.l 
l2009h . 

In this paper we address the issue of hemispherical 
asymmetry by estimating the power spectrum in the two 
hemispheres by using the QML: our application of QML in 
this context is novel and extremely useful since the aggres- 
sive masking needed to reduce residual foreground contami- 
nation might be e ven more problematic for polarization than 
for t emperature (|Bunn et alj 120031 ; ISmith and Zaldarriagal 
l2009h . 

Our main aim is to test whether other asymmetries in 
full temperature-polarization pattern are present around the 
most recently determined axis defi ned by the direction (0 = 
107, 4> = 226) (|Hansen et alJl200Sh . where 9 and are the 
Galactic colatitude and longitude, respectively. 

This paper is organized as follows. In Section 2 we de- 
scribe our methodology by reviewing the algebra of the QML 
estimator. We also discuss the data set used and introduce 
the R and D estimators, the ratio and the difference of the 
power in the two hemispheres respectively. In Section 3 we 
discuss our results including the related Monte Carlo un- 
certenties based on 1000 simulations. We discuss Planck 
predicted performances in probing the hemispherical asym- 
metries in Section 5, while in Section 6 we draw our main 
conclusions. 



2 DESCRIPTION OF THE ANALYSIS 

2.1 Angular Power Spectra Estimation 

In order to evaluate the angular power spectra, we use the 
BolPol code, a QML estimat or. The QML formalism was 
introduced in iTegmarkl dl997h and exten ded to polarization 
in lTegmark and de Oliveira-Costal (|200lh ■ In this section we 
describe the essence of the me thod. Further details can be 
found in lGruppuso et ail (|2009T ) where BolPol) has been ap- 
plied to WMAP 5 year low resolution data. 

Given a map in temperature and polarization 
x = (T,Q,U), the QML provides estimates Cf - with X 
being one of TT, EE, TE, BB, TB, EB - of the angular 
power spectrum as: 
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where the Fisher matrix Fi y , is defined as 
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(signal plus noise contributiorQ) and Cf is a fiducial power 
spectrum. 

Although an initial assumption for a fiducial power 
spectrum Cf is needed, the QML method provides unbi- 
ased estimates of the power spectrum contained in the map 
regardless of the initial guess, 



(Cf) = Cf . 



(4) 



where the average is taken over the ensemble of realizations 
(or, in a practical test, over Monte Carlo realizations ex- 
tracted from Cf). 

On the other hand, the covariance matrix associated to 
the estimates, 



(ACfACf') = (F- 1 )^*' . 



(5) 



with C = S(Cf) + N being the global covariance matrix 



does depend on the assumption for the fiducial power spec- 
trum Cf: the closer the guess to the true power spectrum 
is, the closer are the error bars to minimum variance. Ac- 
cording to the Cramer-Rao inequality, Eq. <(Sj tells us that 
the QML has the smallest error bars. We thus call the QML 
an 'optimal' estimator. 



2.2 Data set and Simulations 

In this Section we describe the data set that we have consid- 
ered and the corresponding simulations we have produced to 
analyze it. We use the temperature ILC map smoothed at 9.8 
degrees and reconstructed at HealPi30 (|Gorski et a.1.1 [2005h 
resolution N a ^ E = 16, the foreground cleaned low resolution 
maps and the noise covariance matrix in (Q,U) publicly 
available at the LAMBDA websit^S We have added to the 
temperature map a ran dom noise realization with variance 
of 1/iK 2 , as suggested in iDunklev etal] (|2009h . Consistently, 
the noise covariance matrix for TT is taken to be diagonal 
with variance equal to lfiK 2 . 

To perform the analysis, we have built the masks for the 
two hemispheres def ined by the direction (6 — 107, <j> = 226) 
l|Hansen et al.l [2009) and combined them with the Galac- 
tic WMAP 5yr low resolution temperature and polarization 
mask. Maps and covariances for the two sky regions (namely 
North and South) have been consistently tailored to the pro- 
duced masks (see Figure (TJ. 

To assess the significance of the power asymmetries 
found in the data, our results have been tested against 
Monte Carlo simulations. A set of 1000 CMB+noise sky 
realizations has been generated: the signal was generated 
from the WMAP 5 year best fit model, the noise through a 
Cholesky decomposition of the global (T, Q, U) noise covari- 
ance matrix. We then computed the angular power spectra 
for each of the 1000 simulations using BolPol and built two 
figures of merit as explained in the next subsection. 



1 Note that, in principle it is possible to include in this matrix 
residuals from foreground subtraction. This is the case for the 
WMAP foreground reduced covariance matrix we employ here- 
after 

2 http://healpix.jpl.nasa.gov/ 
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Table 1. Probabilities (in percentage) to obtain a smaller value 
than WMAP low resolution data for TT angular power spectrum 
and the D-estimator. 

D A£ = 2-8 A£ = 2-16 At = 2-32 Al = 2-40 



TT 



86.2 



96.9 



99.8 



99.1 






Figure 4. TT. Number counts (y-axis) versus R (x-axis) for the 
range Al = [2,40], The vertical line stands for the WMAP 5 yr 
data. The probability to obtain a smaller value than the WMAP 
one is 99.8%. 



Figure 1. Mollweide projection of the observed Northern (red) 
and Southern (blue) hemisphere at iVgide = 16. The (light blue) 
circle in the Southern hemisphere corresponds to the region of 
the Cold Spot, whereas the light grey region corresponds to 
the WMAP low resolution galactic mask for temperature (upper 
panel) and polarization (lower panel). 



2.3 Estimators 



We define the following quantities 
x 1 



C 



N/S 



(In 



1) 



E 



£(i + 1) f,X,N/S 



(6) 



where Cf' N and Cf' s are the estimated angular power spec- 
tra obtained with BolPol observing only the Northern (W) 
and the Southern ('5") hemisphere respectively, outside the 
galactic plane. As above, X runs over the spectral types. 

Two estim ators can be built as follows: the ratio R x , 
as performed in lEriksen et all |2004h . 



TyX f-iX l /~iX 
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and the difference D J 
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of the two aforementioned quantities. In the following, we 
will drop the index X for R and D specifying only the spec- 
trum we refer to. 

For our application to WMAP data, both estimators 
have been considered for TT, while only the D estimator 
has been applied to the other spectra (EE, TE, BB, TB 
and EB), because of unfavorable signal-to-noise ratio of the 
WMAP data in polarization. 



3 RESULTS 

The six angular power spectra TT, EE, TE, BB, TB, EB 
are presented in Fig. [2] and [3] Our results for TT, shown in 
th e upper panel of Fig. [ 21 are consistent with those obtained 
bv lEriksen et all l|2004 ). 

In Fig. 2] we show the R estimator distribution for the 
range £ — 2 — 40. For this estimator we obtain that the 
probability of having the WMAP value is as lo w as .2%, 
which agrees with the results by lEriksen et al. | (|2004T) . In 
Table[T]the probability of obtaining the WMAP value for the 
D estimator is computed for the following four multipoles 
ranges: 2 - 8, 2 - 16, 2 - 32 and 2 - 40. See Fig. 5 for the full 
empirical (Monte Carlo) probability distribution functions. 
Note that the R and D estimator detect a comparable level 
of anomaly in the multipole range 2 — 40. 

In Table [5]we provide results for polarization and cross- 
spectra. As mentioned above, only D is considered and com- 
puted for the four aforementioned multipoles range. The es- 
timator R, in fact, is not well defined any time the denom- 
inator approaches to zero, which might be the case for 
highly noisy spectra. Although Table [5] does not show any 
significant deviation from the symmetry for polarization and 
cross-spectra, it is nonetheless worth noting the behaviour 
of EE in the range £ = 2 — 16, for which the probability 
of having the WMAP value is as low as 3.5%, and of BB 
in the range £ = 2 — 8 where the probability decreases to 
2.2%. Moreover, an unexpected statistics seems to show up 
for TE in the range 2 — 40 where the probability of having 
the WMAP value is 0.5%. However, this mainly comes from 
the multipoles between 32 and 40, which are close to the 
threshold of reliability of the QML on N B id a = 16 maps. 

We also report on the possible contribution to the 
North-South power asymmetry given by the Cold Spot 
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Figure 2. QML estimates for TT (upper panel), EE (middle panel) and TE (lower panel) from WMAP 5 year N s ^ e = 16 maps. Solid 
(blue) line is for the angular power spectrum of the Southern hemisphere (blue region of Fig. [TJ, whereas dotted (red) line is for the 
Northern one (red region of Fig. [TJ- Dashed line shows the WMAP 5 year best fit, taken as fiducial power spectrum for the analysis. For 
reference, we also show the error bars of the QML computed from a Monte Carlo of 1000 sky realizations of the Northern hemisphere with 
the global (T,Q,U) noise covariance matrix (error bars from the Monte Carlo on the Southern hemisphere are basically undistinguishablc 
from the ones plotted). 




Figure 3. As in Fig. [2] but for BB, TB, EB (from top to bottom) 
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Table 2. Probabilities (in percentage) to obtain a smaller value 
than WMAP low resolution data 



D 


Al = 2-8 


At = 2-16 


At = 2-32 


Al = 2-40 


TE 


59.9 


16.9 


75.6 


99.5 


EE 


5.4 


3.5 


28.8 


34.3 


BB 


97.8 


79.5 


71.9 


81.0 


TB 


80.9 


54.1 


42.8 


91.7 


EB 


61.3 


54.6 


74.4 


21.8 




Figure 5. TT. All the panels present the number counts (y-axis) 
versus the D statistics (a;-axis), being the latter in units of fj,K 2 . 
For each panel, D has been computed within a different range of 
multipoles. Top left: Al = [2,8]. Top right: Al = [2, 16]. Bottom 
left: At = [2, 32]. Bottom right: Al = [2, 40]. Vertical lines for the 
WMAP 5 yr data. 



found by IVielva et all (|2004h (see also ICruz et al l |2005l . 
l2007h ). By masking out the Cold Spot (see light blue spot of 
Fig.[T]) with a circle of radius 8 degrees - a conservative choice 
compared to its size of 5 degrees - we have not found any 
significant deviation from the Ct obtained without masking 
it out. This might be due to the fact that the low resolution 
of our data set prevents us from exploring properly the an- 
gular scales of interest for the physical size of the Cold Spot 
(I ~ 40). Moreover, the smoothing process applied to the 
temperature map might be also responsible for washing out 
features like the Cold Spot. Nonetheless, a possible connec- 
tion between these large scale anomalies has been claimed 




Figure 6. As in Fig.[5]but for EE. 






1 







Figure 7. Comparison between WMAP and Planck for EE 
(upper panels) and for TE (lower panels). Left panels: number 
counts for D with Al = [2, 16]. Right panels: D with Al = [2, 8]. 
In all the panels, dark blue is for WMAP and light blue is for 
Planck. Units for D in the x-axis are [iK 2 . 



recently in Bernui (2009), where a different estimator with 
respect to the one adopted in the present analysis has been 
exploited. Withstanding all the caveats set forth above, our 
analysis suggests that the Cold Spot has little to do with 
WMAP 5 year asymmetries. 



4 PLANCK FORECASTS 

The Planck satellite (|The Planck Collaboration! l2006t ) has 
been launched on May 14th, 2009, and it will measure CMB 
anisotropics with unprecedented precision. In order to assess 
its capabilities in probing the hemispherical asymmetry, we 
consider the nominal sensitivity of the Planck 143 GHz 
channel, taken as representative of the results which can 
be obtained after the foreground cleaning from various fre- 
quency channels. The 143 GHZ channel has an angular res- 
olution of 7.1' (FWHM) and an average sensitivity of 6fj,K 
(HAfiK) per pixel - a square whose side is the FWHM size 
of the be am - in temperature (polarizatio n), after 2 full sky 
surveys (|The Planck Collaboration! 120061 ). We assume uni- 
form uncorrelated instrumental noise and we build the cor- 
responding diagonal covariance matrix for temperature and 
polarization, from which, through Cholesky decomposition 
we are able to extract noise realizations. 

As expected, we notice that no significant improvement 
will be achieved with Planck for the TT spectrum, since 
both WMAP and Planck are cosmic variance limited for 
the range of multipoles considered here. On the contrary, 
polarization and cross-spectra do benefit from the Planck 
increased sensitivity. In Fig. 13 we plotted our forecasted 
distribution for the estimator D for TE and EE on top of 
the same distribution for the WMAP case: the shrinking of 
the distribution due to the higher sensitivity is more than 
evident. Moreover, for these two cases we find that also the 
R estimator yields valuable information. Finally, we do not 
expect to be able to apply the R estimator for any spectrum 
involving B because of the low level of signal, if any. Dis- 
tributions of D for BB, TB and EB are analogous to what 
shown for TE and EE. 



6 F. Pad, A. Gruppuso, F. Finelli, P. Cabella, A. De Rosa, N. Mandolesi, P. Natoli 



5 DISCUSSIONS AND CONCLUSIONS 

Using an optimal power spectrum estimator, we have con- 
firmed the power asymmetry for TT found by Eriksen et al.l 
(12004 ) along the direction reported in lHansen et al, I l|2009h . 
see Table [1] and Fig. 5. Considering the same axis, we have 
extended such analysis to the other spectra (TE, EE, BB, 
TB and EB) considering only the estimator D, defined in 
Eq. l|8jl. because the noise level of WMAP permits the use 
of R (see Eq. ©) only for TT. 

Since our implementation of the QML (|Gruppuso et al.l 

2009) is capable of handling the full noise covariance matrix 
in (T,Q,U), the analysis of the present paper is joint for tem- 
perature and polarization. The information encoded in CMB 
polarization is complementary to the temperature and is im- 
portant to t est for possible asymmet r ies in polarization (see 
for instance iDvorkin. Peiris and Hul l|2008l ) for the descrip- 
tion of the polarization field in models that break statistical 
isotropy locally through a modulation field). We confirm the 
TT anomalies that have been already reported by severals 
groups. Our analysis of polarized and cross-spectra does not 
show significant anomalies, as from Table [2] 

The origin of the these hemispherical asymmetries is 
still unknown. They can be primordial or due to so me resid- 
ual fo reground or systematic effect. For instance, in lLi et all 
(2009) an anomalous correlation between temperature and 
observation number has been claimed to be present in the 
WMAP 5yr data, potentially impacting the large scale pat- 
tern of CMB maps (including the power asymmetries). In 
that paper, this effect has been related to an imbalance in 
the differential observation scheme of the WMAP mission. 
Planck is observing the sky with a totally different scheme 
and therefore is free from this particular systematic effect. 

Planck will be able to furtherly confirm the tempera- 
ture anomalies and shed new ligh t onto the polarization sec- 
tor: t he quality of Plan c k data dThe Planck Collaboration! 
20061: iBersanelli et all l20ld : iMandolesi et all l20ld : 
Mennella et al.l 120101 ) is expected to be sufficiently high to 



use the R estimator even for polarization and cross spectra. 
We show the improvements for the D estimator expected 
from Planck in Fig. [7] 

No significant differences have been found by masking 
the Cold Spot in the Southern hemisphere with a disk of 8 
degrees of radius (conservative choice). Despite of the possi- 
ble causes we mentioned in Section 3, we think that our re- 
sults in this respect is worthy of note. Further investigation 
is needed since a correlation between these two anomalies 
(i.e. the Cold Spot an d the North- South asymmetry) has 
been recently claimed (|Bernuil [20091 ). especially in the po- 
larization se ctor where the pro perties of the Cold Spot are 
still unclear (jVielva et alj|20ich . 
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